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Specific heat and /i^^SR measurements in Gd(hfac)3NITiPr molecular magnetic chains: 
indications for a chiral phase without long range helical order 
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Low temperature specific heat C(T) and zero- field muon spin resonance (/x+SR) measurements 
were performed in Gd(fifac)3NITiPr, a quasi one-dimensional molecular magnet with competing 
nearest neighbor and next-nearest neighbor intrachain exchange interactions. The specific heat data 
exhibit a A-peak at To=2.08±0.01K that disappears upon the application of a 5 Tesla magnetic field. 
Conversely, the /i+SR data do not present any anomaly at T « 2 K, proving the lack of divergence 
of the two-spin correlation function as required for usual three-dimensional (3D) long range helical 
order. Moreover, no muon spin precession can be evinced from the /x+SR asimmetry curves, thus 
excluding the presence of a long range ordered magnetic lattice. These results provide indications 
for a low T phase where chiral order is established in absence of long range helical order. 

PACS numbers: 64.60.Fr, 75.40.Cx, 76.75.+i, 75.50.Xx 



I. INTRODUCTION 

For three-dimensional (3D) XY helimagnets, 
Kawamur^ proposed a new universality class as a 
consequence of frustration. In addition to the SO (2) 
symmetry of the spin variable Sr, there is a further 
degree of freedom, the spin chirality,^ defined as 
Ki2 = [Sri X Srj]^/! sin((5a)| (Q is the pitch of the 
helical ground state, and Sra are spins on nearest 
neighbor planes perpendicular to Q). For T = 0, the 
order parameter K12 = ±1 describes the clockwise or 
anti-clockwise degeneracy Z2 of the helical structure. If 
chiral order and spin order occur simultaneously, the 
order parameter space becomes Z2X S0(2), leading to a 
new universality class. 

For quasi one-dimensional (ID) helimagnets, a dif- 
ferent behavior was predicted by Villain."^ In addition 
to the low temperature (T < Tjv) 3D long range heli- 
cal phase and to the high temperature (T > Tq) para- 
magnetic one, Villain proposed that a chiral phase is 
stable for Tn < T < Tq. The chirality order pa- 
rameter associated to a nearest neighbor spin pair is 
K„ = [S„ X Sn+i]"^/! sin((5a)|. The chiral phase can 
be described as a collection of corkscrews that turn 
all clockwise or all anti-clockwise, whereas their phases 
are random;"^ Due to the Ising nature of chirality, Tq 
results greater than T^v because the chirality-chirality 
correlation function is stronger than the usual spin- 
spin correlation. In fact the chiral correlation length 
defined via the chirality-chirality correlation func- 



tion (kik„+i) ~ Aexp(— ria/^s), diverges exponentially 
as temperature is decreased ~ exp(|J|/T) (where 
J is the exchange constant), while the spin correlation 
length ^s, defined via the spin-spin correlation function 
(Si • S„+i) « A' exp(— na/^s), diverges as a power law 
« \ J\/T.^ As a consequence, an anomaly is expected at 
To for the physical properties related to (ki«;„+i), such 
as the magnetic specific heat, while the quantities re- 
lated to (Si • S„+i), such as the magnetic susceptibility, 
are expected to behave critically only at Tjv. At Tq, a 
continuous phase transition or a first-order one can be 
present,- 

Villain's prediction'^ was tested in the molecular- 
based quasi-lD magnet Gd(hfac)3NITiPr^'^ [hfac is 
hexafluoro-acetylacetonate and NITiPr is 2-isoPropyl- 
4,4, 5 ,5-tetramethy 1-4 , 5-dihydro- 1 H-imidazolyl- 1-oxyl 3- 
oxide]; iPr are isoPropyl organic radicals (with spin 
s = 1/2) that alternate with rare-earth Gd'^^ magnetic 
ions (with spin S — 7/2) along the chain direction, z. 
The spins interact through competing nearest-neighbor 
( Ji > 0) and next nearest-neighbor ( J2 < and < 0) 
exchange constants along the chain 
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J2(S ) + £'(5'f 2„-l) — 5M-BH • Si^2n-1 



(1) 



where the index n refers to the spin position along the 
«-th chain. Z? > is an effective anisotropy favouring the 
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spins to lie in the xy plane (perpendicular to the chain 
axis, z). H is an external magnetic field and g,g' are the 
gyromagnetic factors of the spins S and s, respectively. 
The next-nearest neighbor antiferromagnetic exchange is 
dominant, leading to a helical ground state with pitch 

±Qa ~ ±cos^-^ 2{S+S') ' where S = {\J2\S^)/{JisS), 

5' = (|J^|s2)/(JisS') and 2{5 + 5') > 1. Different chains 
interact through a weak ferromagnetic interchain ex- 
change (0 < Jj_ << Ji) 
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so that the total Hamiltonian of a collection of 
M weakly interacting chains is 7i = Winter + 
^iLii'^i-intra)- Magnetic susceptibility and zero field 
specific heat measurements^ were systematically per- 
formed on Gd(hfac)3NITiPr: while no anomaly was 
found in the magnetic susceptibility, a A-type peak, like 
a second order phase transition, was clearly observed 
in the specific heat for T « 2K. Moreover, a linear 
T-dependence was found for C(T) at very low tem- 
peratures, inconsistently with helical 3D long range or- 
der (which would give T^). This overall behavior was 
interpreted^ in terms of a chiral phase transition for 
To w 2K although the phase transition at Tn to he- 
lical long range order, predicted by Villain,^ was not 
observed down to the lowest investigated temperature 
(Tinin = 0.175K for Gd(hfac)3NITiPr). This was tenta- 
tively explained^ taking into account aging effects due to 
the moderate lability of the organic radicals: even in very 
good and fresh samples, chains are reduced to segments 
of finite length and this has the effect of destroying the 
phase coherence necessary for the onset of 3D helical long 
range order. On the contrary, chiral order is compatible 
with the presence of finite segments since, for its onset, it 
is only required that the chirality of the segments is the 
same, a much less strict condition to be fulfilled. In fact, 
while the A-type anomaly in the specific heat is very in- 
tense and robust and its position does not change even in 
degradated samples, aging effects strongly influence the 
magnetic susceptibility>S 

In this paper, we report new specific heat measure- 
ments in magnetic field that demonstrate the magnetic 
nature of the phase transition at T « 2K. In order to ob- 
tain further information about the behavior of the two- 
spin correlation function, muon spin resonance (/i"'"SR) 
experiments in zero-field were also performed. Such a 
technique is particularly useful since it does not require 
the application of an external magnetic field, which in 
a sample with helical long range order is expected to 
cause the onset of complicated spin arrangements ( "fan" 
phase) J To our knowledge, no /i'^SR work on insulating 
hclimagnets has been reported in the literature to date. 

The paper is organized as follows. In Section II the 
new specific heat data are reported, while in Section III 
the results of //"''SR experiments are shown. Finally, the 
conclusions are drawn in Section IV. 



FIG. 1; The specific heat (normalized to the gas constant 
R) versus temperature of Gd(hfac)3NITiPr, measured in zero 
field (filled circles), in a field of 0.5T (T=Tesla)(open circles), 
and 5T (open squares). 



II. SPECIFIC HEAT MEASUREMENTS 

Details about the preparation of samples can be found 
in Ref . 8. The heat capacity of different pellets of pressed 
Gd(hfac)3NITiPr micro-crystals was measured by using 
three different calorimeters and techniques: the adia- 
batic, the ac— and the relaxation method. The heat sig- 
nal was minimized in order to reduce the temperature (T) 
variation with respect to the equilibrium state. The best 
results, on which we concentrate in the following, were 
obtained by using the relaxation method for which the T- 
modulation was less than 0.5% of the base temperature. 
This implies an experimental resolution of logj^g ^ ~ ^2.5 
in terms of the reduced temperature t = \T/To — l\. Mea- 
surements in magnetic field were performed by using a 
Quantum Design PPM-7T System with a •^He insert and 
we begin presenting the new results of these experiments. 

Fig. 1 shows the specific heat C, normalized to the gas 
constant R = 8.314J mol^^ K~^, as a function of tem- 
perature. In zero field a A- anomaly is clearly visible at 
2.08±0.01K while, at lower temperatures, C/R scales to 
zero with a quasi linear T-law. Note that the value of 
the transition temperature is reproducibly found in dif- 
ferent samples (compare also data in Ref. 0|) in spite of 
their different quality and of the height of the A-anomaly 
. The latter gets smoother in a magnetic field of 0.5T 
(T=Tesla) and is completely removed by the applica- 
tion of a 5T magnetic field. This definitively proves the 
magnetic origin of the transition, a fact that could not 
be surely deduced from results of previous experimentsi^ 
Concerning the order of the phase transition, an impor- 
tant feature is the absence of latent heat: we evaluated 
an upper limit of 7xl0~^ J mol~^. This fact seems to 
indicate that the phase transition is of the second order 
in zero field, although the roundness of the peak, which 
can be due to the polycrystalline nature of the sample, 
prevents from determining the value of the critical expo- 
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nent aiS It is worth noticing that the rounded peak at 2K 
and the absence of a latent heat might be justified even 
in the case of a first-order chiral transition by the effect 
of microscopic random quenched impuritiesiSiii that are 
surely present in the magnetic chain. 

In order to gain information about the magnetic order, 
we determined the magnetic entropy removed in the low 
temperature region T < Tq, and in particular that below 
the anomaly. (The magnetic entropy for T > To was not 
evaluated since for Gd(hfac)3NITiPr the lattice contribu- 
tion cannot be estimated in a reliable way owing to the 
complexity of the crystal structure and to the absence of 
a non-magnetic isomorph.) Such a calculation was per- 
formed using the specific heat data relative to sample CI 
in Fig. 2 of Ref. C, since it is considered the best sample 
hereto obtained, owing to its sensibly higher values of 
C/R. We found that the residual magnetic entropy be- 
low Tq is rather large: 59% of the total magnetic entropy, 
which is given by R[ln{2s + 1) + ln(2S' + 1)] = 4i?ln2pi^ 
Such a large value of the residual magnetic entropy be- 
low To is probably due to a high concentration of de- 
fects that prevents the onset of a strong short range or- 
der in the paramegnetic phase, as it usually happens in 
quasi one-dimensional systems. For the sake of compar- 
ison, we remind that in the quasi one-dimensional anti- 
ferromagnet TMMC, the residual magnetic entropy be- 
low the Neel transition temperature T/v was estimated 
to be nearly 1% of the total magnetic entropy>i^ To ob- 
tain the entropy removed below the anomaly at Tq, one 
has to subtract the spin wave contribution (responsible 
for the linear T-dependence of C observed well below Tq) 
from the total specific heat. In this way one obtains the 
value AS — 0.59R for the magnetic entropy related to 
the anomaly, in fair agreement with the value expected 
for an Ising system ASising = -Rln2 « 0.69i?. Such a 
result may be interpreted as a further indication for the 
onset, at Tq = 2.08K, of a chiral transition which removes 
the Z2 degeneracy between clockwise and anticlockwise 
chirality. 



III. fi+SR. MEASUREMENTS 

The /L("'"SR data were collected in pellets of pressed 
Gd(hfac)3NITiPr micro-crystals at the ISIS facihty 
(Rutherford Appleton Laboratory, UK), in the tempera- 
ture range 0.3-300 K, in zero magnetic field (for data in 
longitudinal applied field, see Ref. 113) . After a proper 
background (titanium sample holder) subtraction, the 
muon asymmetry A (whose total value is A^0.21, close to 
0.23, the expected value on the EMU beamline at RAL), 
can be analyzed by means of a single "stretched" expo- 
nential decay e"*^*/'^^ in the whole T-range (see Fig. 2). 

This fact implies, in the whole investigated T-range : 

(i) a distribution of relaxation rates A, witnessed by 
the stretched behavior; 

(ii) the absence of a local field due to long-range mag- 
netic order at the muon site (in presence of long-range 



FIG. 2: Muon asymmetry at two different temperatures. The 
solid line is a fit of the T—4K data (open squares), obtained 
by means of a stretched exponential. As can be evinced from 
the figure, the stretched behavior is also followed by the data 
at r=0.3K (open circles). 

magnetic order, one should observsi^ an oscillation in the 
asymmetry A, caused by the precession of the muon spin 
around the local field itself, or a loss of asymmetry) ; 

(iii) as in the case of short-range magnetic order 
the muon asimmetry should display a Kubo-Toyabe-like 
behavior, one guesses a wide distribution of local fields 
values at the muon site, further broadened by powders 
average. 

Further support for the absence of long-range order 
comes from the analysis of the temperature behavior of 
the muon longitudinal relaxation rate A (see Fig. 3). 

The parameter A represents the muon relaxation rate 
as a result of the muon-lattice interaction. The general 
formula for A in a weak collision approach for fiuctuations 
faster than the muon resonance frequency, lor, isi^ 

Acx72[J+(c^fl) + J„(c^fl)] (3) 

where 

J±= J e-*"«*(h±(t).h±(0)) (4) 

are the spectral densities, at the resonance frequency lor, 
of the correlation functions for the local hyperfine (dipo- 
lar and contact) field components h±. Let us now assume 
that the electronic spin components Sa (coming from the 
Gd ion) and Sa (coming from the iPr organic radical) 
have isotropic fiuctuations. 

To analyze the muon relaxation rate behavior in terms 
of scaling laws, we consider the three types of correlation 
functions that involve two spins. ^ Thus, Eq. Q can be 
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FIG. 3: Muon longitudinal relaxation rate measured in 
Gd(lifac)3NITiPr as a function of temperature in zero applied 
magnetic field. 

rewritten in terms of the collective spin components Sq 
and Sq as 

= E Elhqn^^qll^.ql^ (5) 

lj = l,2 q 1 

where Iji — S and (72 = s; Fq is the decay rate of the 
collective spin fluctuations, while hq is the Fourier trans- 
form of the lattice functions that couple the muons to the 
spins of the magnetic ions. Now, we will show that, in 
the event that long range order sets in below a transition 
temperature Ttr, then a divergence in the muon longitu- 
dinal relaxation rate A is expected (at least for H — 0) 
upon approaching Ttr from higher temperatures. In fact, 
by expanding the hq vector around the value correspond- 
ing to the critical wave vector Q (which characterizes the 
long range order below the transition temperature Ttr), 
and by using conventional scaling arguments for the q- 
depcndence of the correlated fluctuations ^i^. Eq. takes 
the form 

where h^g are the static fields at the muon sites, result- 
ing from the ordered configuration of the Gd and radical 
spins, while qi are the wavevectors measured with respect 



to the critical ones; £_i are the correlation lengths; are 
the Heisenberg exchange frequencies. The decay rate of 
two-spin fluctuations is given by Tq. = {t^e/£,i )g{qi£.i), 
where g{qiS,i) and f{qi^i) are homogeneous functions of 
the product x — qiS,i. By transforming the q-summation 
into an integral, and noting that the latter converge to a 
number of the order of unity, one finally obtains 

By setting cx {T — Tti)~'^\ a divergence of the form 

3 

A cx ^ (T - Ttr)— (8) 

2=1 

with rii = Vi{zi — 1 — 77i), is expected. As can be eas- 
ily seen from Fig. 3, the muon longitudinal relaxation 
rate A does not present any divergence at T « 2 K for 
Gd(hfac)3NITiPr, thus indicating that this temperature 
does not mark any transition Ttr to a long range order 
for the two-spin correlation function. 

Finally we observe that one of the main motivations 
of the present /^"""SR measurements was that of investi- 
gating the behavior of the muon longitudinal relaxation 
rate A around the temperature for which the specific heat 
presents a peak and that the main result of the /i"'"SR in- 
vestigation presented in this paper is the absence of any 
anomaly in A at Tq ~ 2 K. The strong increase of A 
observed below 1 K may have different possible explana- 
tions: e.g., a progressive increase of the coherence lengths 
^i's or a resonance with a spin- wave mode. Further mea- 
surements of A for temperatures lower than IK should 
be necessary in order to ascertain the origin of such a 
feature, but this goes beyond the scope of the present 
work. 



IV. CONCLUSIONS 

In conclusion, we have performed specific heat and 
/i+SR experiments in the molecular magnetic chain 
Gd(hfac)3NITiPr. While the four-spin correlation 
function (probed by the magnetic specific heat) showed 
a seemingly continuous magnetic phase transition at 
T K, 2K, the two-spin correlation function (probed by 
/i^SR) did not show any anomaly related to a transition 
to conventional long range helical order (similar to 
that present in the 3D helimagnets Tb, Yb, Dy, for 
instance).^ Moreover, below T « 2K no precession of 
the muon spin or loss of muon asymmetry similar to 
what happens in the case of long-range magnetic order, 
was observed. Taking also into account that at very low 
temperatures (0.175-0.8 K) a linear T-dependence of the 
magnetic specific heat was observed in Gd(hfac)3NITiPr 
(a clear indication for isotropic ID magnetic behavior) 
and relying upon an estimation of the residual magnetic 
entropy below the anomaly, we are led to conclude that. 
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below To = 2.08 K, an ordered chiral phase without 3D 
long range helical order sets in. 
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